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NONLINEAR  INTERNAL  FLOW  ANALYSIS 

The  aim  of  this  report  is  to  report  the  findings  in 
nonlinear  internal  flow  convection  heat  transfer  analysis.  One  of 
the  major  types  of  loads  which  a  high  speed  aerospace  structure 
encounters  is  due  to  high  thermal  fluxes  especially  on  the 
leading  edges  and  the  associated  thermal  structural  interactions. 
Often  these  type  of  fluxes  are  associated  with  nonlinear 
phenomena  like  compressible  effects,  viscous  effects,  turbulence, 
etc.  and  the  only  possible  solution  to  these  is  through  numerical 
methods.  The  finite  element,  finite  difference,  and  finite  volume 
methods,  and  to  some  extent  boundary  element  method  are  among  the 
most  popular  numerical  methods  which  have  found  widespread  use  in 
computational  mechanics.  Each  of  these  methods  has  its  own 
advantages  and  disadvantages.  The  finite  element  methods  is 
mathematically  very  elegant,  can  handle  complex  geometries 
easily,  and  can  model  very  complex  boundary  conditions.  The 
finite  difference  and  finite  volume  methods  on  the  other  hand, 
are  conceptually  very  simple  and  straight  forward  to  implement  on 
computers .  The  advent  of  powerful  computers  and  work  stations  has 
resulted  in  these  numerical  methods  becoming  very  powerful  design 
tools.  Despite  all  their  advantages,  the  results  of  these  methods 
are  very  heavily  dependent  on  the  uncertainty  of  the  input  data. 
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For,  after  all,  numerical  methods  are  approximate  solutions  to 
differential  equations  (both  ordinary  and  partial) ,  which  in  turn 
are  often  approximations  to  physical  models.  Hence  the  quality  of 
the  numerical  solution  depends  strongly  not  only  on  how  good  a 
particular  numerical  method  is,  but  also  on  how  accurate  the 
physical  model  is.  A  very  obvious  example  of  this  is  numerical 
solutions  of  turbulence  equations.  Because  of  their  stiff 
nonlinearities,  one  has  to  often  resort  to  several  tricks,  the 
most  extreme,  yet  popular  one  being  the  adjustment  of  the  values 
of  turbulence  constants..  This  is  often  required  to  get  a 
converged  solution,  the  accuracy  of  which  is  not  known  to  the 
designer. 

Space  reentry  vehicles  like  ICBM's,  manned  space  crafts,  and 
space  shuttles  are  subjected  to  excessive  heat  fluxes  at  the  time 
of  space  reentry.  One  method  of  preventing  them  from  reaching 
hign  temperatures  would  be  to  let  the  outer  surfaces  which  are 
subjected  to  such  extreme  fluxes  ablate.  While  this  solves  the 
problem  of  overheating  and  protects  the  internal  structure  of 
these  vehicles,  it  also  changes  the  shape  of  these  vehicles  and 
makes  them  unsuitable  from  an  aerodynamic  viewpoint.  This 
disadvantage  can  be  overcome  by  using  internal  ablators  protected 
by  a  fixed  shape  outer  radiating  structure.  During  periods  of 
excessive  heating  internal  ablation  occurs  and  the  gaseous 
products  transpire  through  the  outer  porous  radiating  structure. 
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In  the  following  pages  a  simple  one  dimensional  analysis  of  the 
thermal  ablation  problem  is  presented.  This  has  the  inherent 
disadvantage  that  the  analysis  is  applicable  at  the  point  on  the 
surface  where  heat  flux  is  likely  to  be  the  highest,  like  the 
stagnation  point.  In  fact  it  is  this  point  on  the  surface  where  a 
one  dimensional  approximation  is  most  suitable.  A  schematic  of 
the  stagnation  point  geometry  is  shown  in  Fig.  T.  It  consists  of 
a  fixed  shape  outer  structurethe  temperature  of  which  is  limited 
to  T  .  An  ablator  is  placed  behind  this  outer  structure.  The 
ablator  is  assumed  to  act  like  a  sink  both  in  terms  of  mass 
supply  and  thermal  energy  supply  so  that  it  can  be  assumed  to  be 
at  a  constant  temperature  Tb.  When  the  surface  temperature  of  the 
ablator  Ts  rises  due  to  heat  radiation  and  reaches  a  value  of  Ta, 
it  ablates  and  gaseous  products  transpire  through  the  fixed  outer 
surface  thereby  reducing  the  heat  transfer  to  the  outer  surface. 
The  heat  transfer  mechanism  between  the  outer  structure  and  the 
ablation  surface  is  by  radiation  across  the  gap.  Of  the  heating 
rate  g0(t)  incident  on  the  outer  structure,  a  part  qr/Q(t)  is 
radiated  outward  while  the  rest  qr/i(t)  is  radiated  inward  to  the 
ablation  surface.  When  ablation  occurs,  the  surface  temperature 
of  the  ablator  increases  to  a  fixed  value  of  Ta. 

GOVERNING  EQUATIONS: 

At  the  porous  shield  the  following  equation  is  obtained  by  a 


3 


thermal  energy  balance : 


q0(t)  =  qr,0(t)  +  qr,i(t) 


(1) 


In  the  above  equation  q0(t)  is  the  aerothermodynamic  heating  on 
the  vehicle  surface,  while  qr,o(t),  qr,i(t)  are  the  radiative 
exchanges  between  the  shield  and  outer  atmosphere  and  inner 
ablator  surface  respectively.  Since  the  outside  atmosphere  can  be 
assumed  to  be  at  absolute  zero,  qr,o(t)  can  be  approximated  as 


qrA*)  =  ev/ffUt) 


(2) 


Qrfi)  = 


TOO  -  tM 

[1  +  fe  -  ew] 


(3) 


where  it  has  been  assumed  in  Eq.  (3)  that  the  heat  transfer 
mechanism  between  the  porous  shield  and  ablator  surface  is  due  to 
radiation.  During  ablation  as  the  gaseous  products  move  through 
the  porous  surface  they  reduce  the  heat  flux  q0(t)  .  By  making  a 
heat  balance  at  the  ablator  surface  the  following  equation  is 
obtained. 


QnetV)  =  qoit )  -  rjrhH0(t )  (4) 

where  H0(t)  is  the  total  free  stream  enthalpy,  and  rj  is  an 
effiencY  for  the  vaporization  process.  During  ablation  at  steady 


k 


state,  a  thermal  energy  flux  babalce  equation  at  the  shield  wall 
is  given  by 

?»(')  -  VhH0(t)  =  ShoTJ«)  +  £  (5) 

e"a  £w] 


While  ablation  takes  place,  Ts(t)  =  Ta.  At  the  ablator  surface 
an  energy  balance  gives  the  following  equation: 


rao  -  T]m 

[1  +  £t  ~  ew] 


=  mHv  +  6bcbj-t[(Ts  -  Tb)6 ] 


(6) 


In  the  above  equation  Hv  consists  of  two  parts,  namely  the 
ehthalpy  of  vaporization  also  called  the  enthalpy  of  ablation, 
and  the  enthalpy  rise  before  vaporization.  For  a  planar  slab, 
integral  thermal  thickness  can  be  used  to  approximate  the  thermal 


conduction  into  the  interior  as  follows: 


"  *(f> =  -  w 


(7) 


Hence  Eq.  (6)  can  be  rewritten  as 


.  ffK(')-^(0]  .  ,  (Ts  —  Tb) 

“Ti  +  -  £j  =  +  XA—W^\ 


(8) 


Thus  for  a  combination  of  geometry  of  the  trajectory,  sheild  and 
ablator  properties  the  behaviour  of  the  system  can  be  obtained  by 
solving  Eqs.  (5),  (6),  and  (8) .  When  no  ablation  takes  place,  the 
mass  ablation  rate  can  be  set  to  zero  and  the  equations  solved 
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for  the  variables  Tw(t)  ,  Ts(t),  and  0(t).  When  ablation  takes 
place,  the  surface  temperature  Ts(t)  is  set  equal  to  a  constant 
Ta  and  now  the  variables  to  be  solved  forare  Tw(t)  ,  m(t)  ,  and 
0(t)  .  The  next  step  in  the  process  ig  to  recast  the  equations  in 
nondimensional  form  which  is  shown  below. 


Non  dimensional  variables  cab  be 

easily  obtained 

normalizing  each 

variable  against  some 

reference  state 

follows : 

II 

Ih? 

Tw,o  ^  Tw  <  TWjmax 

II 

^1^ 

Tb  - 
■t<ts<  1 

1  a 

&  ^  ^max 

(9) 

Hv 

'<=  mqmtm 

ft  ^  ft  max 

II 

rh 

z<r, 

and  F(r)  =  ^ 

m  <  i 

In  the  above  set  of  equations  qm  is  the  maximum  input  heat  flus 
of  the  trajectory  which  occurs  when  t  =  tm.  Using  these  non 
dimensional  quantities,  Eqs  (5),  (6),  and  (8)  can  be  reformulated 

in  dimensionaless  form.  With  no  ablation,  we  can  set  the  mass 
flow  rate  term  to  be  zero  and  write  equation  (5)  as 
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(10) 


F{x)  =  RwTw(r )  -  RaTs(r) 
where  Rw  ans  Ra  are  the  radiation  parameters  defined  as  follows: 

(11a) 

016) 

with 

q*r  =£w  a  r*4 


Rw  =  1^(1  + — — - ) 

Qm  £a  £a£w 


R  -  (l£\4  n  ,  £g 

qm  T*>  +  Ea  +  Ew-  EaEJ 


We  can  use  Eq.  (10)  to  solve  for  the  temperature  ratio  to  get 
rt(r)  =  ^[F(r)  +  Rafs(r)]  (12) 

Using  the  above  in  the  non-dimensional  form  for  Eqs .  (6),  and  (8) 

we  get 

F(T)  +  (£a-w  -  1)  Rafs( T)  =  ~  Tb)d]  (13) 

ea,w  F(z)  +  (ea>w  -  1)  Rafl(r)  =  -  Tb)]  (14) 

with 


o  Qj£bTa@g 

Pa  qmtm 


* 'b£a 

qrrfi  g 


(15) 
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(16) 


£a,w  — 


2£q  4*  £w  £<j£v 


Equations  (12)  -  (14)  are  a  set  of  three  nonlinear  coupled  set  of 
equations  which  need  to  be  solved  for  Ts,  Tw,  and  6.  A  further 
simplification  can  be  made  to  the  set  of  equations  (10),  (13), 

and  (14)  if  a  temperature  function  i|>  is  introduced  as  follows: 
m  =  I7»  -  Tb]2  (17) 


and  an  expression  for  the  derivative  for  ij; 

d  _  BJW  ±  A  fT 
dry’  AQf!>\ 


can  be  written  as 
(18) 


where  B [op]  ,  C[ip],  and  A  are  defined  as 
Bhp]  =  eatW  F(t)  +  (ea>w  -  1)  Ra(Tb  +  ftp)* 

Cbp]  =  Ea,w  F( x)  +  (sa>w  -  1)  Ra(Tb  +  ftp)\Tb  -  ftp) 

A  =  fia/ a 


(19) 

(20) 

(21) 


and  F'  =  dF/dx.  A  numerical  integration  of  Eq.  (18)  gives  us  the 
funciton  ip(x)  .  Once  the  temperature  function  \p(x)  is  known,  we 
can  easily  get  Ts(t),  the  nondimens ional  surface  temperature  from 
Eq.  (17).  Using  Ts  (x)  ,  we  van  get  Tw  from  Eq.  (12)  and  S(x)  from 
Eq.  (14)  respectively. 


On  the  other  hand  when  vaporization  takes  place  as  when  the 
surface  temperature  of  the  ablator  reaches  the  ablation 
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temperature  Ta,  mass  flux  must  be  taken  into  account.  Thus  we 
have  from  Eq.  (5)  after  nondimensionalizing. 


?t(T)  =  4-[^(r)  ~  rj/XHoir)  +  Ra] 

where 


(22) 


HJt)  =  Ml 

Since  the  material  is  ablating,  we  have  Ts  =  1 .  The  dimensionless 
forms  of  eqs .  (6)  and  (8)  are 


'«,»  Hr)  +  -  1)  R„  =  (1  +  +  /3„(1  -  Tb 

ea,»  Hr)  +  -  1)  R„  =  +  ea,^H0)^  +  y„° 

By  rearranging  Eq.  (24)  we  get 

dp  _  gg,w  F(r)  +  (ea,w  -  1)  Ra  -  ^(1  -  Tb) 

dT  (lt  +  ea,wvHo) 


(23) 

(24) 


(25) 


A  single  differential  equation  for  5  is  obtained  by  combining 
Eqs.  (23)  and  (24) 

dd  _  gg,w  F(t)  +  (£a,w  ~  1)  Rg(  1  ~  IT)  ,Ya,n 

dr  Ml  ~  Tb)  W<5  (2( 

where  parameter  IT  is  defined  as 


Yj  _  14"  £  a,w1jH0 

1  ~  Ahv  .  Yf 

ft  4"  £a,w^]H0 


(27) 


It  has  been  pointed  out  by  Camber os  (1989,  1996)  that  for  typical 
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Q-folation  materials  II  ~  1.  Equation  (26)  is  solved,  numerically 
for  6  and  the  information  is  used  in  Eq.  (25)  to  solve  for 
dimensionless  mass  flux.  Presently  calculations  are  underway  to 
verify  the  one  dimensional  equations.  For  the  case  of  two 
dimensional  problem,  the  resulting  set  of  equations  would  be 
integral  equations . 

NOMENCLATURE : 


A  =  dimensionless  ablation  parameter 

c  =  specific  heat  at  constant  pressure,  J/kg-K 

F  =  dimensionless  heating  rate  function 
Ahv  =  enthalpy  of  ablation  J/kg 

q  =  heating  rate  per  unit  area,  W/m2 

R  =  dimensionless  radiation  parameter 

T  =  temperature  K 

t  =  time  s 

V  =  velocity  m/s 

x,y  =  coordinates  m 


Greek  Symbols 

(3  =  dimensionless  parameter 
Y  =  dimensionless  parameter 
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=  dimensionless  thermal  thickness 
e  =  emissivity 

Q  =  density  kg/m3 

K  =  thermal  conductivity  W/m-K 

H  =  dimensionless  mass  flux 

a  =  Stephan-Boltzmann  constant  =  5 . 67  x  10-8  W/m2K4 

x  =  dimensionless  time 

=  dimensionless  temperature  function 

SUBSCRIPTS : 

a  =  ablator  property 

b  =  ablation  material 

o  =  stagnation  point  value 

r  =  radiative 

s  =  ablation  surface 

w  =  porous  shield  wall  value 
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Fig.  1:  One  dimensional  geometry  near  stagnation  point  of  ablator  surface 
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